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INTRODUCTION 


One  type  of  specimen  commonly  used  in  tests  designed  for  the 
determination  of  a  material's  fracture  toughness  can  be  looked  upon  as  a  double 
cantilever  beam,  as  shown  in  Figures  1  and  2. 


Figure  2.  Compact  fracture  toughness  specimen. 
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The  data  collected  la  such  tests  to  determine  the  material's  fracture 
toughness  property  Includes  the  gap  opening,  2d',  at  a  given  load,  P' ,  and  for 
a  given  crack  length,  V  .  The  stress  distribution  In  such  specimens  Is  a 
complex  one  and  so  Is  the  resultant  strain  distribution  and  its  contribution  to 
the  gap  opening.  The  computation  of  that  portion  of  the  gap  opening  attributed 
to  the  elastic  bending  of  a  cantilever  beam  is  the  subject  of  this  report. 

These  calculations  are  based  on  an  elastic  beam  theory.  The  shortcomings  of 
the  computed  gap  opening  derived  are  as  follows: 

1.  The  fracture  rarely  propagates  without  attaining  a  plastically 
deformed  region  near  the  tip  of  the  crack. 

2.  The  calculations  do  not  account  for  the  deformation  in  the  remaining 
portion  of  the  specimen,  namely  at  the  areas  x'  >  V  ,  beyond  the  crack  length. 

The  added  gap  opening  due  to  the  propagation  of  plastic  deformation  should 
be  computed  separately  using  the  technique  used  by  this  author  elsewhere  (ref 
1).  The  contribution  of  the  deformation  in  the  region  x’  >  V  towards  the  gap 
opening,  can  also  be  computed  separately.  (Neither  one  of  these  calculations 
is  Included  here.) 

In  treating  each  half  of  the  specimen  as  a  tapered  cantilever,  one  has  to 
realize  that  Its  neutral  axis  in  bending  is  not  parallel  to  the  specimen's 
plane  of  symmetry  (x*  axis  In  Figure  1),  but  rather  at  an  angle  a  to  it  (as 
will  be  shown  later,  the  effective  neutral  plane  deviates  from  the  bisecting 
plane  of  each  cantilever,  too).  As  a  consequence,  the  computation  of  the 

^Boaz  Avitzur,  "Retained  Deflection  in  Circular  and  Concentrically  Hollowed 
Beams  After  Local  Removal,"  to  be  published. 
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beading  moment  sad  of  the  shear  forces,  as  well  as  that  of  the  resultant 
stresses  thereof,  will  be  based  on  the  component  P  -  P' •  cos  a  of  the  applied 

h© 

load  P',  where  P  acts  normal  to  the  axis  of  symmetry,  x  (where  x  -  [x*  - -  tan 

a]  cos  a). 

PARALLEL  STRESS  COMPONENTS 

As  a  first  approximation,  the  gap  opening  6  ■  $bend  +  ^shear  is  calculated 
in  the  cartesian  coordinates  system  x-y-z,  which  is  at  an  angle  o  to  the 
x’-y'-z*  coordinates  and  only  the  normal  component,  Oxx»  atu*  shear 
component,  o^  ■  o^,  are  considered.  This  is  followed  later  by  correcting  the 
stress  components  to  be  parallel  to  the  external  surfaces  at  the  boundaries 

h0  +  2x*tan  a 


Figure  3.  Definition  of  cartesian  coordinates  and  major  parameters. 
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The  basic  dimensions  in  the  x-y  coordinates  and  their  relation  to  those  in 
the  x'-y'  coordinates  (see  Figure  3)  are  as  follows: 

V 

h  m  -  =  the  beam' s  cross-sectional  height  at  the  point  of  load  application 

cos  a 


^o 

x  «  (x'  -  —  tan  a)*cos  a,  where  x'  is  being  measured  along  the  top  surface  of 
2 

the  beam,  which  is  the  fracture  toughness  specimen,  along  the  crack 
surface. 


I 

t 

P 

P 

h 


ho* 

*  (!'  -  — -  tan  a) ’cos  a  =  distance  of  crack  tip  from  point  of  load 
application 

*  2* tan  a 

»  P' *cos  a  =  component  of  load  normal  to  plane  of  symmetry 

*  P'*sin  a  =  component  of  load  parallel  to  plane  of  symmetry 

■  h0  +  t*x  =  beam's  height  normal  to  the  plane  of  symmetry  at  a  distance  x 
from  the  point  of  load  application. 

Thus,  the  cross-sectional  moment  of  inertia  becomes 


bh3  (hQ+tx) 3 

12  12 


(1) 


and  the  local  curvature  p  of  the  beam  is 


1  M(x)  Px  P  * 

P  E*I(x)  E’I(x)  Eb  (hQ+tx)3 

and  the  total  deflection  due  to  bending  moment  is 
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7o 
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0odx  +12 


Li*]* _ I _ 

Eb  0  0  <h0  +  3 


dx’dx 
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(3) 
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for  which 

P  ,x  x  PI 

6  -  eo  +  12  —  J - 5  *  dx  -  0O  +  12  - =  [ - 

0  Eb  o  (hQ+tx) 3  0  Ebt2  hQ  + 


P  1  “O  x 

Ebt2  ^  hQ  +  tx  2(hQ+tx)2^o 


80  -  12  —  * - -5 

0  Eb  2h0(h0+tx) 2 


where  0  *  angle  of  deflection 

0Q  *  angle  of  deflection  at  the  point  of  load  application. 
However,  at  the  tip  of  the  crack,  x  =■  l,  0£  »  0.  Therefore, 


0p  -  0  =  Qn  -  12  — 


Eb  2h0(h0+t i)  ■ 


from  which  we  can  determine  that 


0O  -  6 


Eb  h0(h0+ti): 


By  applying  Eq.  (5)  to  Eq.  (3),  one  gets 


y  m  Vo  +  F  eo<»*  +  12  ~  J*  J*  """-I  dxdx 
o  Eb  o  o  (hQ+tx) J 


P  x  x  xz 

y0  +  6 - {- - 5  f  dx  -  /  - ■=  dx} 

Ebh0  (hQ+ti)2  o  o  (h0+tx) 2 


P  l2  x  1  ho  x 

y°  +  6  {7Z~1Z~o2^  1  ~  :i  tho+tx  "  2ho  log(hQ+tx)  -  - - - ]  } 

Ebh0  (h0+t2/)  o  t3  hQ  +  tx  o 


P  l2  1  2  ho  +  tx  hQ  +  tx 

y° +  6  sr  {”7r2 x  -  r*  [r""r  * tx  ■  2h° 108  — 1 }  (6) 


Ebh0  (h0+ti)‘ 
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0,  one  gets 


Since  at  the  tip  of  the  crack,  x  -  SL,  and  the  deflection  y  - 

SL 


o  m  yQ  +  6 - { 


Ebh0  (hQ+t  i) 2  t3  ^  +  tl 


1  2ho+ti 

Ir— — “  tt  -  2ho  log  — - - ]  }  (7) 


hQ  +  ti 


However,  for  a  straight  cantilever  where  t  *  0,  one  has  to  determine 


7l‘° 


2h0  +  tx 


hn  +  tx 


y  +  6 - { - -z  -  lim(— ;  [ - tx  -  2h0  log 

Ebho  (h0+ti)2  t-KD  fc3  ho  +  tx 


-])> 


where 


1  2h0  +  tx  hQ  +  tx 

lira  {~t  [ - tx  -  2h0  log - ]  } 

t+0  t3  hc  +  tx  hQ 

from  which 


xJ 

lim - -s 

t-K)  3(h0+tx)  1 


(8) 


p  o3  p3  p  o  3 

lim  y0  =*  -6 - {—• - - -}  =  -4  --  (--)  =  yQ' 

t>0  Ebh0  h02  3h02  Eb  hQ 

which  is  the  same  as  being  reported  in  handbooks  (ref  2)  for  a  straight 

cantilever  of  a  uniform  cross-section. 

Another  method  of  computing  the  gap  opening  is  by  invoking  Castigliano 
theorem  (ref  3),  equating  the  internal  deformation  energy  with  that  of  the 
outside  work  done  on  the  structural  member.  This  method  will  facilitate 
computation  of  the  gap  opening  due  to  bending  moment  and  due  to  shear  force 
(and  any  other  forces  that  Impose  internal  deformation  work). 

The  distribution  of  that  component  of  normal  stresses,  oxx>m,  which  is 
parallel  to  the  axis  of  symmetry  and  which  results  from  the  bending  moment, 
as  follows: 


(9) 


’s 


is 


^Raymond  J.  Roark  and  Warren  C.  Young,  Formulas  for  Stress  and  Strain,  Fifth 
Edition,  McGraw-Hill,  NY,  1975,  p.  98. 

3A.  C.  Ugural  and  S.  K.  Fens ter,  Advanced  Strength  and  Applied  Elasticity, 
American  Elsevier,  NY,  1981,  pp.  146-148. 
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M  P 

axx,m  -  -  *  y  -  -  '  x-y  (10) 

where  P  »  P*  cos  a,  a  *  arctaa  (t/2)  ■  arctan  (Ah/2 Ax),  and  I  *  bh3/12 

■  “  *  (V-tx)3. 

The  component  of  shear  stress  which  lies  in  a  plane  normal  to  the  plane 
of  symmetry  and  in  the  direction  of  the  axis  of  symmetry  (or  opposite  to  it), 
°yx  m  and  *ts  equivalent  0Xy>ra»  is  assumed  to  be  due  to  gradient  in  that 
component  of  the  normal  stress  doxx  m/dx,  resulting  from  the  bending  moment. 
This  is  computed  as  follows: 

h/2  b/2  Aaxx,m 

2  /  /  — - - dzdy 

o  o  Ax 

rT  m - ■ - - • - 

2 

where  M 

°xx,m  =  ~  ‘  y 

Thus 

daxx,m  d  M  d  12  P  d  x 

dx  dx  I  dx  b(h0+tx) 3  b  dx  (hQ+tx) 3 

where 

d  x  ^o  ”  3tx  hg  —  2tx  i  tx 

dx  (hQ+tx)3  (hQ+tx)4  (hQ+tx)4  (hQ+tx)3^  hQ  +  tx^ 

and  thus 

P  (h0+tx)/2  b/2  h0  -  2tx 

24  -  j  /  - -  y*dz*dy 

k  y  °  (hQ+tx)  i2P  (h0+tx)/2  hD  -  2tx 

°yx,m  - - - - - J  7:t:~74  *  ydy 

b  by  (hD+tx) 4 

3P  hQ  -  2tx  3P  ho  “  2tx  y 

°vx  m  ”  —  * - Z  [ ( hn+tx)  2  -  4 y 2 ]  -  —  • - -  [1  -  (2 - )  21 

’  2b  (h0+tx)4  2b  (h0+tx) 2  hQ+tx 

3P.  hQ  -  2tx  y 

.  - — •- - -  [1  -  (2 - )  2]  *  cos  a  (11) 

2b  (hQ+tx) 2  h0+tx 
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Thus 


”-3  -  +  ?  «■>  C  ~p  «. 

where 

x  x2  i  hQ  +  tx  2h0  +  3tx 

^ o  (h0+tx)3  dX  1 3  ^  8  h(}  2(h0+tx) 2 

Thus 

12  9  1  P2  hQ  +  tx  2h0  +  3tx 

ubend3  “  <“3  +  5  ">  ^  [1°8  “  ^h^txj'3  tXl  “  P<5bend3 

from  which 

12  9  1  P  hQ  +  tx  2h0  +  3tx 

«be„d3  -  (-3  +  -  ;)  ~  -----  -  ---“--2  tx]  (26) 

and  as  before 

12  9  1  P  +  tx  2h0  +  3tx 

«'bend3  -  +  ;  ->  5  U*  ~~~  -  £~~p  tx]  -cob  a  (26a) 

At  the  limit,  as  t  +  0,  for  a  beam  of  a  constant  cross-section, 

P  12  9  1  ho  +  t3c  2ho  +  3tx 

llm  ® bend 3  *  ”  lim  +  “  ->U-°8  ~ - o7u~Z~~7~2  txl  } 

t+0  bE  t+0  tJ  5  t  h0  2(h0+tx) 2 

where 

d  hQ  +  tx  2h0  +  3tx 

—  [log - -  tx] 

l  hQ  +  tx  2h0  +  3tx  dt  ho  2(h0+tx)  *■ 

lim  {--  [log - -  tx]  }  -  llm - 

t+0  t  hQ  2(h0+tx) 2  t+0  d 

db  (t3> 


x  [(2h0+3tx)x  +  3tx2]<hD+tx) 2  -  2(2h0+3tx) *tx*(h0+tx) *x 

hQ  +  tx  2(h0+tx)4 

■  lim - - - 

t+0  3t2 


21 


Since  <Jxxjm  (or  <xrr>a)  is  tensile  in  the  range  0  <  8  <  a  (or  0  <  y  < 
(h0+tx)/2)  and  compressive  in  the  range  -a  <  0  <  0  (or  -(hQ+tx)/2  <  y  <  0), 
ayy  m  m  °xx,m*can  9  maintains  its  upwards  component  and  these  do  not  cancel 
each  other  (Figure  5).  Nevertheless,  no  discontinuity  prevails  at  0  *  0  (or 
y  *  0) ,  since  Oyy>m  *  0  on  that  (internal)  surface.  The  contribution  of  the 
stress  component,  Oyy>m  to  the  deflection  is  accounted  for  by  the  energy 
balance  in  its  totality. 

Deflection  Due  to  Bending  Forces  Only 
Since 


^xx.m 


P  xy 

12 - -5 

b  (hQ+tx) 3 


Jrr,m 


12 


p  J  t^2  *y 

b  (hQ+tx)2  (hQ+tx) 3 


from  which  the  local  strain  energy,  due  to  bending  moment,  ubencj3  becomes 


rr,m 


ubend3  *  °rr,m'err,m  *  “  *  “2^  + 


tV 


E  b^  (hD+tx) 2  (hD+tx)6 


Hence,  one  can  compute  the  total  strain  energy  due  to  the  bending  moment  as 
follows: 


P2  <x  (h0+tx)/2  ,b/2 
ubend3  *  J  Ubend3’<lv  “  576  — -  /  J 


xV2 


xV2 


b*E  00 


,  [1  + - -5]  - --dz-dy-dx 

o  (hQ+tx) 2  (hQ+tx) 6 


P2  ,x  ,(hc+tx)/2  t2y2  x2y2 


where 

(h0+tx)/2 


288  rr  /  J 

bE  o  o 


t2y2  x^2  x2  y3  t2  y9  (hg+tx)/2 


[1  + - - — -5) - - — -*dy*dx 

(hQ+tx) 2  (hQ+tx)6 


[1  +  (h0+tx)zl  (h0+txj6  dy  ”  (h”+tx>9  lr +  (h0+tx>2  r'o 


.2  .,3 

6 


JL  1  c 

( hQ+tx) ^  l24  +  1601 
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of  any  cross-section  at  a  distance  x  from  the  line  of  load  application  is 

L  -  L'+x  -  ~  +  x  from  the  origin  of  these  cylindrical  coordinates  (where 
t 

t  -  2 ‘tan  a),  and  the  distance  l  of  any  point  (x,y)  on  a  plane  x  and  at  a 


Figure  5.  Radial  and  cartesian  components  of  normal  stresses. 
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CONVERGING  STRESS  COMPONENTS 

By  assuming  that  the  normal  stresses,  ffn,  are  parallel  to  the  plane  of 
symmetry,  one  violates  the  boundary  condition  at  y  ■  ±  (h0+tx)/2,  in  a  tapered 
beam*  Therefore,  it  will  be  assumed  here  that  the  normal  stresses,  arr>  at  the 
surfaces  of  y  ■  ±  (h0+tx)/2  (or  9  -  ±  a)  (see  Figure  4)  are  parallel  to  these 
surfaces  and  that  the  radial  stresses  gradually  align  themselves  with  the  axis 
of  symmetry  (8-0).  Or,  in  other  words,  these  are  radial  stresses,  as  their 
subscript  suggests,  in  a  cylindrical  coordinate  with  the  origin  at  L'  ** 
(h0/2)*tan  a  -  h„/t  distance  from  the  point  where  x  -  0.  Thus,  the  distance  L 


Figure  4.  Conversion  from  cartesian  to  cylindrical  coordinate  systems. 


.*•  /-  >  . 


.•  v  V  *  •  * 
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fc.* 


Hence  ,  , 

1+v  F  12  1  ho  +  2ho  +  5t*  t3  ho  +  tx  .... 

6shear2  E  b  1  5  t  14  8  h(J  2(h0+tx)2  24  ^ 

and  , 

1+v  P'  12  1  h0  +  tx  2ho  +  5tx  t3  N,  +  tx 

6,  ,  „  til  _.f4.log  -- - 3  — - - r  tx]  +  —  log - >cos2  a 

shear2  E  b  5  t  hQ  2(h0+tx) 2  24  ho 

(24a) 

Combining  Eqs.  (22)  and  (24),  one  gets 

1+v  P  12  1  hQ  +  tx  2h0  +  5tx  t3  hQ  +  tx 

- { - [4 -log - 3 - -5  tx]  +  -7  log  — 7 - > 

«3h«r2  _  E  b  5t _ _ X*o+‘*> _ “ _ _ 

5bend2  P  12  ho  +  tx  2ho  +  3tx  t  hD  +  tx 

—  {-7  [log  — . — -j  +  7  log  ““7 - 1 

bE  t3  hQ  2(h0+tx) 2  4  hQ  (25) 

As  before  (in  computing  lim  <$bend) » 

t>0 

t3  ho  +  tx 

lim  {—  log - )  -  0 

t+0  24  ho 

Therefore 

12  P  1+v  x 

ltm  68hear2  •  “  *  7  (”") 
t+0  5  b  E  ho 

as  In  Eq.  (16),  and 

^shear2  3  h© 

iim - (l+v)(— )  2 

t+0  ^bend2  5  x 

as  in  Eq.  (17). 


17 


^LisLtJLtSL  •*■  +■ 


f.  > 


Thus,  the  total  shear  strain  energy  due  to  the  combined  bend  and 
compressive  forces  and  over  the  entire  volume  of  the  beam  Ushear2  becomes 
(h0+tx)/2 


^shear2  *  2b  0 


o  o 


1+  v  P2  .x  (h0+tx)/2  „ 

ushear2*^y’dx  9(h0  24h0txf4t2x2) ' 

Ci  D  O  O 


1  8y2  +  16y^  ^ 

^(h0+tx)4  (hQ+tx)^  (h0+tx) ® 


t4  1  4 | y |  4y2 

+  —  f - t - ;  + - r]  }dydx 

4  (h0+tx) 2  (hQ+tx) 3  (h0+tx) 4 


As  In  the  derivation  of  Eq.  (4) 


(h0+tx)/2 


I 


8y; 


I6y 


^ (hQ+tx) 4  (h0+tx)^  (h0+tx)®*  7  15  (hQ+tx)^ 


and 

(h0+tx)/2 


1  4 1  y  |  4y  1111  1 

^(hQ+tx)2  (hQ+tx)3  (h0+tx)4^  7  ^ 2  2  6^  Iiq  +  tx  6(h0+tx) 


Thus 


1+v  P2  x  12  N>2  "  4hotx  +  4t2*2  t4  1 


l+v  P*-  x 

ushear2  -  ----  -  J  {• 


E  b 


o  5 


(h0+tx) 


24  hQ  +  tx 


■}dx 


As  in  the  derivation  of  Eq.  (14) 

^  hQ2  -  4h0tx  +  4t2x2  1  h0  +  tx  2ho  +  5tx 


(h0+tx) 


-  dx  -  -  (4 -log - 3 - „tx ) 


2(h0+tx) 


whereas 


fx  dx  __  1  1  ^  +  tx 

o  ho  +  tx  t  hg 


Thus 

1+v  P2  12  l  h0  +  tx  2h0  +  5tx 

^s hear 2  -  ~  “  {-•-•[4-log  —  - 3  — — - 

2  b  5  t  hQ  2(h0+tx) 


ho  +  tx 

2  tx]  +  —‘log  — ; - ) 


f 

24 


ho 


PS, 


shear2 
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as  in  Eq.  (9),  since 


t  h©  +  tX , 

lia  {-  log - }  ■  0 

t-K)  * 


By  adding  the  two  components  of  shear  stresses  due  to  the  bend  forces,  P' 
cos  a,  and  due  to  the  compressive  forces,  P*  sin  a,  from  Eqs.  (11)  and  (18), 
respectively,  one  gets 


P 

—  {3 


hQ  -  2tx 


y  2  fc2  lyl 

•y*  ‘  °r*-*  T  “»*•!> '  S'  “  (h„+t«) 2  [1 ' <2  1  '  io^Sj'11  ‘  2 


Thus,  the  local  combined  shear  strain  energy,  ushear,  due  to  the  bending  and 
compressive  forces  is 


1+v  p2  —  2tx 


ushear2  *  °yxeyx  “  2  g“  °yx2  “  E“  ^2  {3  (h^)2*[1  "  (2  h^t?  ^ 


c  _  I _  2 

2(h©+tx)  ^  h©  +  tx^ 


1+v  P2  ho2  **  ^h©tx  +  4t2x2 

T  2b2  {9  ~  (h©+tx)  4 


[1-2(2  — )  2  +  (2  — )4] 
h©+tx  h©+tx 


h©  -  2tx 


-  3 - r  [1  -  (2  — 

(h©+tx)^  h©+tx 


h©+tx  4(h©+tx) 2 


[i  -  4 - - - +  4  — - - -J  } 

hQ  +  tx  (h©+tx) 2 

where,  due  to  symmetry  around  the  x-axis,  the  term 

ho  “  2tx  v  I y I 

3  Trm*11  ■ <2  rr")2],[1  * 2 

(h©+tx)J  h©+tx  h©  +  tx 

cancels  Itself  on  both  sides  of  the  axis  of  symmetry. 
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«bend2  "  2b  /!  / 


(ho+tx)/2  p2  ,  (h0+tx)/2  x2y2 

o  ubend2*<ly^  “  C  !o  076  (hQ+tx)- 


(h0+tx) 

where,  as  in  the  derivation  of  Eq.  (12) 


dy*dx 


In  addition, 


(h0+tx)/2  x2y2  ^  J.  _ x2 

o  (h0+tx) 6  ^  *  24  (hQ+tx) 3 


(h0+tx)/2  _ 1 2 _  __  t^_ 

o  (h0+tx) 2  ^  2(h0+tx) 


P2  x  x2  t2 

U.n/iJ  -  L.  r  (24 - 5  + - }dx 

hend2  2bE  JQ  (hQ+tx) 3  2(h0+tx) 


P2  24 


- [log(h0+tx)  + - - - -~2]  +  “  log(h0+tx)  } 

2bE  t3  ^  hQ  4-  tx  2(h0+tx) 2  2  c 


P2  12  hQ  +  tx  2h0  +  3tx  t  ho  +  tx 

-  —  [log  — - - ~ tx]  +  -  log  - - >  -  p&bend2  (21> 

bE  t3  ho  2(h0+tx) 2  4  ho 

from  which 

P  12  ho  +  tx  2ho  +  3tx  t  h0  +  tx 

5bend2  “  —  {”5  [log - 5  tx]  +  T  lo8 - J  (22) 

beadZ  bE  t3  h0  2(h0+tx) 2  4  ^ 

and  as  before 


p»  12  ho  +  tx  2ho  +  3tx 


h„  +  tx 


w  1108  ~Z 


- -  tx]  +  -  log - 

2(h0+tx) 2  4  ho 


•  }co8 2  a  (22a) 


which,  as  was  previously  shown  (following  the  derivation  of  Eq.  (12)) 


lia  6bend2  “  “4  “  3 


Considering  the  symmetry  of  the  bean,  one  arrives  at 


-Pt2  lyl  -P’t  lyl 

o  -  -  [1-2 - ]  -  -  [1“2 - ]  *sin  a  (18a) 

^’P  4b(  he+tx)  ho  +  tx  4b(h0+tx)  N,  +  tx 

By  adding  the  two  components  of  normal  stresses  due  to  the  bend  forces, 

P* ‘cos  a,  and  due  to  the  compressive  forces,  P'  sin  a,  one  gets 


°xx  "  axx,m  +  °xx,p 


^  P  xy  P  t  P  {^4  xy  ^  _ 1 


b  (hQ+tx) 3  2b  hQ  +  tx  2b  (hQ+tx) 


hQ  +  tx 


The  term 

{24 - -  -  t} 

(h0+tx) 2 

in  Eq.  (19)  is  not  symmetric  in  respect  to  the  plane  of  symmetry  (y  -  0) 
(except  for  beams  of  constant  cross-section,  where  t  -  0) .  Thus  the  beam's 
geometrical  plane  of  symmetry  is  not  necessarily  the  plane  of  symmetry  for  the 
stress  distribution. 

As  in  the  derivation  of  Eq.  (12),  the  local  strain  energy  due  to  normal 
stresses,  ubend2,  is 

°xx2  P2  x^2  xy  ,  1 

ubend2  ”  °xxexx  - -  — T~  i576  - -  -  48  - -5  t  +  t*> - 5 

E  4b 2E  (ho+tx)1*  (h0+tx) 2  (hQ+tx)2 


However,  in  Integrating  over  the  entire  volume  and  due  to  the  symmetry 


around  the  x-axls,  the  form 


xy 

48 - .  t 

(h„+tx) 2 


cancels  itself  on  both  sides  of  this  axis.  Thus,  the  total  strain  energy  due 
to  the  normal  stresses  and  over  the  entire  volume,  Hbead,  is 


12  P  1+v  x  12  P  1+v  x 

ss  Wl  ■  1 6  T  <4-3)V  '  1  b  T  %>  <l6 

and  together  with  Eq.  (9)  (where  yQ  »  ^bend) 

12  P  1+v  ^x 

^shearl  5  b  E  hQ  3  h0  m 

Urn - - - -  -  (l+v)(— )  2  (17 

t+0  ^bendl  P  x  3  ^  x 

4  --  (— ) 

bE  h0 

Deflection  Due  to  Combined  Bending  and  Compressive  Forces 

If  one  Incorporates  the  effect  of  the  compressive  stresses 

P'  P  sin  a  Pt 

°xx,p  ”  8  n  “  "  b(h0+tx)  cos  a  2b(h0+tx) 

then  the  change  In  compressive  stresses  due  to  change  In  cross-sectional  area 


becomes 


2b  (h0+tx)2  b  (h0+tx)2 


from  which  the  associated  shear  stresses  will  be 


(h0+tx)/2  k/2  ^axx,p  Pt2  (h0+tx)/2  h/2 


—  *dz  *dy  -2  — - -5  / 

k  2b(h0+tx) 2  y 


fb/2 

J  dz*( 


-Pt2  ho  +  tx  -Pt2 

'  H  ^  Ak/k  4-1 


2b(h0+tx) ‘ 


4b(ho+tx)  ^  ho  +  tx^ 


-P’t.  y 

- —  [1-2 - ]  «sln 

4b(ho+tx)  ho  +  tx 
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and  as  before 

12  P»  1+v  1  ho  +  tx  2ho  +  5tx 

fi,«haarl  - - {4  log - 3  - ;  tx >co# *■  a 

shearl  5  b  E  t  ^  2(h0+tx) 2 

which,  together  with  Eqs.  (12)  and  (12a),  yields: 


^shearl 


12  P  1+v  1  h0  +  tx  2ho  +  5tx 

- —  {4  log - 3  - 5  tx) 

5  b  E  t  ho  2(h0+tx)  shearl 


pi  ho  +  tx  2ho  +  3tx 

12 - -  {log - -5  tx) 

bE  t3  hQ  2(h0+tx) 2 


shearl 


°bendl 


xrr  v  u 

-  {4  log  - 

5  ho 


1+v  h0  +  tx  2ho  +  5tx 

- {4  log - - - 3 - -  tx) 


2(h0+t*>  2  6' shear l 


1  ha  +  tx  2ha  +  3tx 

Uog - -  tx) 

t2  h0  2(h0+tx) 2 


0  bendl 


However,  at  the  limit  as  t+0,  which  Is  equivalent  to  a  straight  cantilever 


of  uniform  cross-section,  one  gets 


hQ  +  tx 


4  log - 

12  P  1+v  ^  2ho  +  5tx 

llm  ^shearl  ‘  ~  '  ~~~  lia  < - - - 3  - - TJ  *> 

t+0  5  b  E  t+o  t  2(hoytx)2 


where 


hQ  +  tx 


+  tx  d  hQ  +  tx 

o  ^  ho 


llm  { - - }  ■  lira  {- 

t+0  t  t+0 


x  x 


d  t+0  h0  +  tx  ho 

r  <*> 

dt 


2h0  +  5tx 


t+0  2(h0+tx) ' 


9P2  4  1+v  p  h0  -  4hotx  +  4t2x2 
ushearl  “  T  15  T  £  (h0+tx)  3  ** 


for  which 


-X  dx  *»o  l  x  N>  1  1  1  2N>  +  tx 

°  o  (h0+tx)3  2t  (h0+tx)2  o  2t  ^ho2  (h0+tx)2^  t  2(h0+tx) 


,x  x*dx 


1  1 


^O  X 


1  4t2x2 


4h  [  41i  w  a  *  -tw  a 

^  Q  (h0+tx)  ^  ^  t*  2  Vi  -VVv  4»Vv\  2^  A  ►  *7  f  Vi  4»^v\  2 


t2  hQ+tx  2(h0+tx) 2  o  t  2(h0+tx) 


,x  xzdx  4 


4t2/  ------  -  -  [log(h0+tx)  + 

o  (h0+tx) 3  t 


hQ+tx  2(h0+tx)  2  o 


4  hQ+tx  2h0  +  3tx 

-  [log - - 

t  hQ  2(h0+tx) 2 


2  txl 


.x  ho  “  4h0tx  +  4t2x2  1  hQ  +  tx 

/  - — - — - ri - • —  dx  -  -  (4  log - + 

o  (hQ+tx) 3  t  ho 


2h0tx  +  t^2  -  4t2x2  -  8hotx  -  12  t*x2 
2(ho+tx)  2 


1  hQ  +  tx  2ho  +  5tx 

-  {4  log - 3 - -  tx} 

t  ho  2(h0+tx) 2 


from  which 


12  P2  1+v  1  h0  +  tx  2ho  +  5tx 

Ushearl  '  “  r  T  ‘  {4  log  - 3  — - -  tx)  -  P 6, 


5  b  E  t 


2(h0+tx) 2 


ahearl 


and  thus 


12  P2  1+v  1  ho  +  tx  2 ho  +  5tx 

6shearl  '  “  ~  7  *  W  log  — - - 3  — - -  tx} 

5  b  E  t  ho  2(h0+tx) 2 


If,  however,  in  the  derivation  of  Eq.  (12)  one  uses  P' *coa  a  for  P  and  is 
to  determine  the  displacement  5'  in  the  x'-f*  coordinate  system,  one  will  get 


the  following: 


hn  +  tx  2h0  +  3tx 


-  12 - ;  (log - o * tx } * coa z  a 

bend  I  bE  t3  6  u  2 


2(h0+tx) 


(12a) 


The  local  strain  energy  due  to  that  component  of  shear  stress  resulting 

from  the  gradient  in  normal  bend  stresses  and  acting  on  planes  normal  to  the 

plane  of  symmetry  and  In  the  axial  direction,  (°r  *t8  equivalent  acting 

on  planes  parallel  to  the  plane  of  symmetry  and  in  the  direction  normal  to  the 

plane  of  symmetry,  Oxy>m)  ushearl,  is 

1+v  „  9P2  1+v  (h0-2tx) 2  y 

ushearl  ’  ^yx.m^yx.m  »  2  ~  0  yx,m  *  2  •  •  T  (h^tx)4"  U  ‘  ^h'+tx*  1 

9P2  1+v  ,  ,  1  8y2  16y4 

2b2  E  °  °  MV(h0+tx)4  (hD+tx) 6  (h0+tx)8J 

for  which  the  total  shear  strain  energy,  Ushear^,  over  the  entire  volume 


becomes: 


9P2  1+v  x  (h0+tx)/2 

ushearl  “  ^2  ~T  *  b*2J  /  (ho^hotx+At^2) 

z  d  is  o  o 


1  8y2  16y4 

^(hQ+tx) 4  (hQ+tx) ^  (h0+tx)3^  7 


for  which 


(h0+tx)/2  i  8y 2  16y4 

o  ^ (h0+tx) 4  (hQ+tx) ^  (ho+tx)0^dy 


y  8y3  16y5  ( 

^(h0+tx)4  3(h0+tx) 6  5(h0+tx)8^o 


16y 5  (h0+tx)/2  4  l 

hQ+tx)  8^o  "  15  (h^+tx) 3 
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Deflection  Due  to  Bending  Forces  Onl 


The  local  strain  energy  due  to  normal  stresses,  %en<n ,  is 


ubendl  “  °xx,m  •  exx,m 


xx, a  P2  x2y2  P2  x^2  _ 

- 144  — s - ?  -  144  —x - 1  *  cos2  o 

E  b2E  (hQ+tx) 6  b^  (h0+tx)6 


Thus,  the  total  bending  strain  energy,  Obeni|i,  over  the  entire  volume  of  the 


(deformed)  beam  becomes: 


x  (ho+tx)/2  p2  v  (1 

ubendl  *  2b/  /  ubend‘drdx  *  288  "  i  / 

o  o  bE  o  o 


P2  *  (h0+tx)/2  x^y2 


(h0+tx) 6 


dy*dx 


where 


(h0+tx)/2  x2y2  1  2  (h0+tx)/2  i  x2 

f  J  <  i 


a  y  a  a 

- -  dy  *  - - t  y  3/ 

(h0+tx) 6  3  (hQ+tx) 6  c 


24  (hQ+tx) 3 


for  which 


IX  -“““3  *  dx  -  "  [log 

o  (h„+tx) 3  t 3 


hrt  +  tx  2h0  +  3tx 


2(h0+tx) 


288  P2  1  ho  +  tx  2ho  +  3tx 

w '  “I  bi  t1 1108  ~z~ '  tT) 


12 - =  {log - 

bE  t3  h„ 


hD  +  tx  2h0  +  3tx 


2(h„+tx) 2 


tx}  ■  P^bendl 


from  which 


PI  hQ  +  tx  2ho  +  3tx 

6bendi  "  12  rr  “3  {log  ~~Z - — 7*2  tx} 

bE  t3  hn  2(h«+tx) 2 


oenai  bE  fc3  *  h<)  2(h0+tx) 2 

which  is  the  same  as  y0  arrived  at  in  Eq.  (7)  through  the  double  integration  of 
the  bending  moment  method,  and  as  such 


P  x 

lim  Sbendl  "  “4  “C") 
t+0  bE  h„ 


as  in  Eq.  (9) 


X 


ZYiq^x.  +  8hotx2  +  St^3  -  4hotx2  -  bt^3 


h0  +  tx 


2(h0+tx)  3 


11m 

t*0 


3t' 


+  4h0tx2  +  It2*3  -  Zho^  -  8h0tx2  -  bt^3  +  4h0tx2  +  bt^3 


lira 

t+0 


3t‘ 

-  Zt2*3^ 

2(h0+tx) 3  a  ^ 

Hi - - - -  lim - -  -  -  (—)  3 

t+0  3t2  t>0  3(h0+tx) 3  3  ^ 


to  which  oae  may  add  that 


1  hQ  +  tx 

lim{-  [log  - 

t+0  t  ho 


2h0  +  3tx 

- --tx]  } 

2(h0+tx) 2 


hQ  +  tx 

log  — - - 

ho  2h0  +  3tx 

lim - 11m - -x  *  x 

t-K)  t  t*0  2(h0+tx)2 


where 


hQ  +  tx 
log - 

h°  x 

lim - *  lim - — 

t>0  t  ^o  +  tx 


x 

ho 


and 

2h0  +  3tx  x 

t+0  2 ( ho+tx) 2  ho 


Thus 

1  h0  +  tx  2h0  +  3tx 

lim  {-  [log - - - - -  tx]  }  ■  0 

t+0  t  ho  2(h0+tx) 2 


Therefore 


6bend3 

t+0 


P  1  x 

12 - (--) 

bE  3  h0 


3 


22 


which,  except  for  the  sign,  is  the  same  as  that  arrived  at  in  Eq.  (9),  and 
which  concurs  with  the  value  given  in  handbooks  (ref  2)  for  a  straight 

cantilever  of  a  uniform  cross-section. 

Even  though  the  shear  stress  approaches  zero  at  the  specimen’s 
boundaries,  y  -  ±  (hQ+tx)/2,  and  thus  it  is  not  in  violation  on  the  boundaries, 
it  is  assumed  here  that,  as  with  the  normal  stresses,  ^yx  the  cartesian 
component  of  a  radial  shear  stress,  OQr,  which  converges  at  (or  radiates  from) 
the  same  cylindrical  coordinate  origin  at  L  *  -(hQ/t  +  x)  as  before. 


/  ty 

r  a  -  /  1  +  ( - ) 

■I11  t, 


h0+tx 


Thus,  since 


!ec,.  -  [i + 


3P  b©  -  2tx  y 

'  2b‘(h^S)2  U  ‘  ‘  1 


3P  h0  -  2tx  y 

*  u  U 


hQ+tx 


Thus,  the  local  shear  strain  energy' due  to  the  bending  forces,  u8t,ear3  is 

1+  v  9p  2  i+  v 

ushear3  “  °0r,me0r,m  “  2  —  <>20r,m  -  ~z  (ho^^tx+At^2) . 


^(h0+tx)4  (hQ+tx) 6  (hQ+tx) 


6  +  /u  ‘ll  +  /u 


(h0+tx) 2 


2 Raymond  J.  Roark  and  Warren  C.  Young,  Formulas  for  Stress  and  Strain,  Fifth 
Edition,  McGraw-Hill,  NY,  1975,  p.  98. 
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from  which  the  total  shear  energy  due  to  bending  forces  only  over  the  entire 
beam's  volume 


^shear3  “  *  r  i 


(h0+tx)/2  ^j,/2 


ushear3  *  dz • dy  *  dx 


9P2  1+v  x  (h0+tx)/2  ,,  1  8y2  16y4 

STlU.  >-1^55  -  ^—“e  *  5^5*1  • 


[1  + - -]  dydx 

(hc+tx)  2J 


where 


(h0+tx)/2  1  8y2  16y4  tV 
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2  24  20  56 
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(hQ+tx)3  l(2  ~  3  +  105  +  (24  20  +  56)t:  1 
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(15  +  1055  (hg+tx) 3 

1  t2  ho  “  4h0tx  +  4t2x2 

-  —  (4  +  — ) - - - 

15  7  (h0+tx) 3 


Thus, 


t2  Pz  1+v  jc  b0  -  4h0tx  +  4t2x2 


5  t4-  Y-  1+v  X 

-  —  (4  +  — ) - r 


- - - dx 


where,  as  in  the  derivation  of  Eq.  (14) 


>2  -  4h0tx  +  4t2x2  1  h0+tx  2h0  +  5tx 

- - - - = - dx  -  -  (4-log - 3  ------“Tl  tx) 
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Hence 


3  4  t  P  1+v  h0  +  tx  2h0  +  5tx 
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shear3  ^  7^  b  E  v>  2fh-+t-v^2 
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and  as  before. 
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Since 


hQ  +  tx  2h0  +  5tx 

lim  (t[4 -log - 3  “--77-72  tx]  }  -  0 

t-+n  hn  2(h_+tx) 


as  in  the  derivation  of  Eq.  (16)  above 

12  P  1+v  x 

lim  5shear3  ”  “  7  ~Z~ 

t-*-0  5  b  E  hQ 

as  it  should  be.  Also, 
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10  t  7 
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2(h0+tx) 2 


h«  +  tx  2h„  +  3tx 


Deflection  Due  to  Combined  Bending  and  Compressive  Forces 


The  same  process  of  considering  the  converging  stresses  used  in  the  above 
calculations  of  the  gap-opening  due  to  the  bending  forces  can  be  used  for  the 
combined  bending  and  compressive  forces* 

As  in  Eq.  (19)  before,  the  Cartesian  component  of  the  combined  forces 


°xx  *  axx,m  +  °xx,p 


P  xy 

—  (24 - 

2b  (hQ+tx) 


ho  +  tx 


Thus,  as  before,  the  full  normal  stress  will  be 


P  /  ty  2  xy  1 

-  ra  =*  —  /  1  +  ( - )  •  {24 - -  -  t  }• - 

’*  2b  hQ+tx  (h0+tx) 2  hQ  +  tx 


and  the  local  strain  energy  due  to  normal  stresses,  ubend4,  will  be 

o2rr  P2  t2*2 


ubend4  ”  °rrerr  ”  ~  *  TTl  l1  + 


E  4b2  (hQ+tx) 2 


,  x2y2  txy  t2 

(576  - t  -  48 - - -  + - -} 

(hQ+tx) 6  (h0+tx) 4  (hQ+tx) 2 


However,  due  to  the  beam's  symmetry,  when  integrating  throughout  the  entire 


beam's  volume,  the  term 


txy 

48 - 

(h0+tx) 4 


from  both  sides  of  the  plane  of  symmetry  cancels  itself,  and  one  gets  for  the 
total  bending  strain  energy,  Ubend4*  the  following: 

x  h0+tx/2  b/2  P2  x.  (ho+tx)/2  t^2 

ubend4  -  /  /  /  ubend4‘dz*dydx  "  TTZ  f  J  t1  +  7r-7r~72]  * 
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where 
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as  has  been  shown  before.  Thus, 
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Since 
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as  before 
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Incorporating  the  bending  and  the  compressive  components  of  the  applied 
force  into  the  computation  of  the  resultant  radial  3hear  stresses,  yields 


the  following: 

p  ho  ■  2tx  y  2  ^  2 
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Thus  the  local  shear  strain  energy  resulting  from  the  radial  shear  stresses  due 
to  the  combined  bending  and  compressive  forces  is 
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However,  when  Integrating  over  the  entire  voluae  one  can  onit  the  tens 
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since  it  cancels  itself  on  both  sides  of  the  plane  of  symmetry*  Therefore 
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and  similarly 
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as  before 
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and 
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SUMMARY 

T  ils  report  provides  equations  for  the  determination  of  that  part  of  the 
gap  opening  in  fracture  toughness  specimens,  which  is  due  to  the  deflection  of 
the  cantilever  portion  of  the  specimen.  The  method  used  is  based  on  beam 
theory  and  it  covers  elastic  deformation  only.  Four  different  field 
combinations  were  considered. 

1.  Where  only  the  stress  components  which  are  parallel  (or  normal)  to  the 
cantilever  beam’s  plane  of  symmetry  are  being  considered. 

2.  Where  the  stresses  are  assumed  to  be  parallel  to  the  beam's  outer 
surfaces,  at  these  surfaces,  and  where  they  gradually  rotate  in  between. 

Each  of  the  above  stress  fields  was  considered  as  the  result  of 

a.  the  bending  component  of  the  applied  force  only,  and 

b.  the  combined  bending  and  compressive  components  of  the  applied  force. 

In  all  of  the  above  cases,  the  contributions  of  bending  stresses  and  shear 

stresses  were  computed.  Table  I  is  a  computer  printout  for  the  computed  half 
gap  opening  for  varying  angles  of  specimen  taper  and  for  varying  crack  length 
to  beam* 8  cross-sectional  height  ratios.  The  gap  openings  to  be  anticipated  by 
each  of  the  above  assumed  stress  fields  are  compared. 
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TABLE  I.  GAP  OPENING  (HALF)  IN  TAPERED  DOUBLE  CANTILEVER  BEAMS 
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TECHNICAL  REPORT  INTERNAL  DISTRIBUTION  LIST 


CHIEF,  DEVELOPMENT  ENGINEERING  BRANCH 
ATTN:  SMCAR-LCB-D 
-DA 
-DP 
-DR 

-DS  (SYSTEMS) 

-DS  (ICAS  GROUP) 

-DC 

CHIEF,  ENGINEERING  SUPPORT  BRANCH 
ATTN:  SMCAR-LCB-S 
-SE 

CHIEF,  RESEARCH  BRANCH 
ATTN:  SMCAR-LCB-R 

-R  (ELLEN  FOGARTY) 

-RA 

-RM 

-RP 

-RT 

TECHNICAL  LIBRARY 

ATTN:  SMCAR-LCB-TL 

TECHNICAL  PUBLICATIONS  &  EDITING  UNIT 
ATTN:  SMCAR-LCB-TL 

DIRECTOR,  OPERATIONS  DIRECTORATE 

DIRECTOR,  PROCUREMENT  DIRECTORATE 

DIRECTOR,  PRODUCT  ASSURANCE  DIRECTORATE 


NO.  OF 
COPIES 


NOTE:  PLEASE  NOTIFY  DIRECTOR,  BENET  WttAPONS  LABORATORY,  ATTN:  SMCAR-LCB-TL 
OF  ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST 


NO.  OF 

NO 

.  OF 

COPIES 

COPIES 

ASST  SEC  OF  THE  ARMY 

COMMANDER 

RESEARCH  &  DEVELOPMENT 

US  ARMY  AMCCOM 

ATTN:  DEP  FOR  SCI  &  TECH 

1 

ATTN:  SMCAR-ESP-L 

1 

THE  PENTAGON 

ROCK  ISLAND,  IL  61299 

WASHINGTON,  D.C.  20315 

COMMANDER 

COMMANDER 

ROCK  ISLAND  ARSENAL 

DEFENSE  TECHNICAL  INFO  CENTER 

ATTN:  SMCRI-ENM  (MAT  SCI  DIV) 

1 

ATTN:  DTIC-DDA 

12 

ROCK  ISLAND,  IL  61299 

CAMERON  STATION 

ALEXANDRIA,  VA  22314 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENG  ACTV 

COMMANDER 

ATTN:  DRXIB-M 

1 

US  ARMY  MAT  DEV  &  READ  COMD 

ROCK  ISLAND,  IL  61299 

ATTN:  DRCDE-SG 

1 

5001  EISENHOWER  AVE 

COMMANDER 

ALEXANDRIA,  VA  22333 

US  ARMY  TANK-AUTMV  R&D  COMD 

1 

ATTN:  TECH  LIB  -  DRSTA-TSL 

COMMANDER 

WARREN,  MI  48090 

ARMAMENT  RES  &  DEV  CTR 

US  ARMY  AMCCOM 

COMMANDER 

ATTN:  SMCAR-LC 

1 

US  ARMY  TANK-AUTMV  COMD 

1 

SMCAR-LCE 

1 

ATTN:  DRSTA-RC 

SMCAR-LCM  (BLDG  321) 

1 

WARREN,  MI  48090 

SMCAR-LCS 

1 

SMCAR-LCU 

1 

COMMANDER 

SMCAR-LCW 

1 

US  MILITARY  ACADEMY 

SMCAR-SCM-0  (PLASTICS  TECH 

1 

ATTN:  CHMN,  MECH  ENGR  DEPT 

1 

EVAL  CTR, 

WEST  POINT,  NY  10996 

BLDG.  351N) 

SMCAR-TSS  (STINFO) 

2 

US  ARMY  MISSILE  COMD 

DOVER,  NJ  07801 

REDSTONE  SCIENTIFIC  INFO  CTR 

2 

ATTN:  DOCUMENTS  SECT,  BLDG.  4484 

DIRECTOR 

REDSTONE  ARSENAL,  AL  35898 

BALLISTICS  RESEARCH  LABORATORY 

1 

ATTN:  AMXBR-TSB-S  (STINFO) 

COMMANDER 

ABERDEEN  PROVING  GROUND,  MD  21005 

US  ARMY  FGN  SCIENCE  &  TECH  CTR 

ATTN:  DRXST-SD 

1 

MATERIEL  SYSTEMS  ANALYSIS  ACTV 

220  7TH  STREET,  N.E. 

ATTN:  DRXSY-MP 

1 

CHARLOTTESVILLE,  VA  22901 

ABERDEEN  PROVING  GROUND,  MD  21005 

NOTE;  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER, 
US  ARMY  AMCCOM,  ATTN:  BENET  WEAPONS  LABORATORY,  SMCAR-LCB-TL, 
WATERVLIET,  NY  12189,  OF  ANY  ADDRESS  CHANGES. 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER, 
OS  ARMY  AMCCOM,  ATTN:  BENET  WEAPONS  LABORATORY,  SMCAR-LCB-TL, 
WATERVLIET,  NY  12189,  OF  ANY  ADDRESS  CHANGES. 
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